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Abstract 

Context. NGC 6397 is one of the most interesting, well observed and theoretically studied globular clusters. The existing wealth 
of observations allows us to study the reliability of the theoretical white dwarf cooling sequences of low metallicity progenitors, to 
determine its age and the percentage of unresolved binaries, and to assess other important characteristics of the cluster, like the slope 
of the initial mass function, or the fraction of white dwarfs with hydrogen deficient atmospheres. 

Aims. We present a population synthesis study of the white dwarf population of NGC 6397. In particular, we study the shape of the 
color-magnitude diagram, and the corresponding magnitude and color distributions. 

Methods. We do this using an up-to-date Monte Carlo code that incorporates the most recent and reliable cooling sequences and an 
accurate modeling of the observational biases. 

Results. We find a good agreement between our theoretical models and the observed data. In particular, we find that this agreement 
is best for those cooling sequences that take into account residual hydrogen burning. This result has important consequences for 
the evolution of progenitor stars during the thermally-pulsing asymptotic giant branch phase, since it implies that appreciable third 
dredge-up in low-mass, low-metallicity progenitors is not expected to occur. Using a standard burst duration of 1.0 Gyr, we obtain that 
the age of the cluster is 12.8 /q 3 j Gyr. Larger ages are also compatible with the observed data, but then unrealistic longer durations of 
the initial burst of star formation are needed to fit the luminosity function. 

Conclusions. We conclude that a correct modeling of the white dwarf population of globular clusters, used in combination with the 
number counts of main sequence stars provides an unique tool to model the properties of globular clusters. 

Key words, stars: white dwarfs — stars: luminosity function, mass function — (Galaxy:) globular clusters: general — (Galaxy:) 
globular clusters: individual (NGC 6397) 


1. Introduction 

White dwarfs are the remnants of the evolution of single stars 
of low and intermediate mass. Consequently, they gather im¬ 
portant information not only about the evolution of their pro¬ 
genitor stars, but also on the properties of their parent popula¬ 
tions. Moreover, since in general terms, the evolution of white 
dwarf stars is relatively well understood, the ensemble prop¬ 
erties of their populations can be used to test the validity of 
phys ical theories that cannot be don e using terrestrial facili¬ 
ties jGarcfa-Berro et alJ 1 1995L1201 lat llsern et aD 1 19921 120081 
l2010t ICorsico et al.ll2012l) . This has been possible because we 
now have a profound knowledge of the physics governing white 
dwarf interiors and envelopes, which has resulted in accurate 
cooling tracks, as well as a wealth of observational data. Even 
more, the quality of the available observational data has al¬ 
lowed us to check the accuracy of our description of the main 
physical proc esses involved in their evolution - see, for in¬ 
stance 1 Althaus et alJ (120101) for a recent review. In this vein, it 
is important to realize that due to the conceptual simplicity of 
the white dwarf populations of globular clusters their proper¬ 
ties can be used to test white dwarf evolutionary models. This 
is particularly interesting because it has been recently shown 
that, under controlled conditions, theoretical white dwarf cool¬ 
ing times obtained using very different numerical codes and 
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techniques do not differ by much. Hence, white dwarf evolution- 
ary times sh ould be at least as reliable as main sequence lifetimes 
(tSalaris et al.ll2013h . 

Several studies have taken advantage of all these prop¬ 
erties, and have determined ac curate ages of both open and 
globular clusters. For instance, [Garcfa-Berro et all d20ld ) de¬ 
termined the age of the old, metal-rich, well-populated, open 
cluster NGC 6791 using t he term ination of the cooling se¬ 
quence (iBedin et al.l 120051. l2008af hl), and solving the long¬ 
standing problem of the discrepancy between the main se¬ 
quence and white dwarf ages. But this is not the only clus¬ 
ter for which we have good age determinations using white 
dwarf evolutionary sequences. Other examples ar e the young 
open clusters M 67 dBellini e t al.ll201 0h. NG C 2158 dBedin et al.l 
I2010I) . and NGC 6819 (BedinefaL 12015 ) , the g lobular clus¬ 
ters M4 ( Hansenet alJ l2002t Bedin et al. 20091). to Centauri 
dBellini et al. 2013 ). and NGC 6397 (Hansen et al.l 201 3h. or 
the metal-rich glob ular cluster 47 Tuc dGoldsburv et al.l 20121 : 
iGarcfa-Berro et aPl2014l) . 

NGC 6397 is the second-nearest globular cluster to the 
Sun, and has been thoroughly observed by the Hubble Space 
Telescope. Thus, for it we have high-quality deep images that 
have allowed us to study not only the lower main sequence, but 
also the white dwarf cooling sequence. NGC 6397 is o ld and 
me tal-poor, being its metallicity [Fe/H]= -1.8 () Han sen et al .l 
1201 3l) . although in the recent lit erature there exist s some dis¬ 
crepancy about its precise value dRicher et al.l[2008 ). The same 
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happe ns for the age of the cluster. For instance, lHansen et ali 
( 2007 1 analyzed the position of the cut-off of the white dwarf lu¬ 
minosity function and, comparing with theoretical coolingjnod- 
els, derived an age T c = 11.47+0.47 Gyr. Similarlv. lWinget et al.l 
d2009t) . fitted simultaneously the main-sequence, the pre-white 
dwarf and the white dwarf regions of the color-magnitude di¬ 
agram, and obtained an age T c = 12.0+jq Gyr. These ages, 
which are based on white dwarf evolutionary models, need to 
be compared with those obtained fitting the l umin osity of the 
main-sequence turn-off (MSTO) - see [Richer et al.l ( 2008) for a 
discussion of different methods and age estimates of NGC 6397. 
In particu lar, it is worth highlighting that using this method 
iGratton et alj _([2003[l derived an age of 13.9 ± 1.1 Gyr. This 
prompted IStrickler et al.l (l2009i) to claim that this age determi¬ 
nation was compatible with the possible existence of a putative 
population of very old helium white dwarfs. On the other hand, 
Anthony-Twaro g& Twarod (l2000h obtained 12.0 ± 0.8 Gyr and 


Chabover et al. (uOOl ) derived 13.4 ± 0.8 Gyr. Thus, the pre¬ 


cise age of NGC 6397 remains unclear, and needs to be in¬ 
dependently evaluated using the most recent white dwarf evo¬ 
lutionary tracks of the appropriate metallicity. To this we add 
that NGC 6397 shares some other interesting characteristics with 
other Galactic globular clusters, such as the unusual lithium en¬ 
hancement of some of their stars - see lPasauini et al.l (12014h and 
referen ces therein - or the possible existence of mul tiple popu¬ 
lations (idi Criscienzo et aPl2010l : lMilone et al.ll2012l) . Thus, de¬ 
riving an independent white dwarf cooling age for this globular 
cluster is of crucial importance. 

Moreover, NGC 6397 is one of the few globular cluster 
for which, in addition to its color-magnitude diagram, we also 
have an accurate white dwarf luminosity function, thus allow¬ 
ing us to inv estigate in detail those issues previously mentioned. 
Specifically, IWinget et al. (12009) used the shape of the white 
dwarf luminosity function of NGC 6397 to constrain important 
properties of white dwarf interiors. They found that to account 
for the observed shape of the white dwarf luminosity function 
of NGC 6397 either a low oxygen mass fraction in the inner 
carbon-oxygen core of typical white dwarfs was needed, or the 
crystallization temperature of the carbon-oxygen dense binary 
plasma should be significantly larger than that of a one com¬ 
ponent plasma. The reason for this is that when the carbon- 
oxygen plasma cry stallizes the oxygen abund a nce in the solid 
phase is enhanced dGarcia-Berro et all 1988a b; Horowitz et al.l 
120101) . This leads to a redistribution of oxygen in the core of 
the white dwarf, because the inner regions are enriched in oxy¬ 
gen, while the outer ones become oxygen-poor. Since oxygen 
is slightly heavier than carbon this chemical separation pro- 
cess releases gravitation a l energy, incr e asing the cooling times 
( Isern et aD 11997112000 1. IWinget et al.l ( 2009 1 found that to fit 
the observations the effects of phase separation upon crystalliza¬ 
tion should be minimized. Hence, either the gravitational energy 
should be released at larger luminosities (or, equivalently, core 
temperatures), or the oxygen abundance in the deep interior of 
typical white dwarfs should be smaller. This, in turn, is an im¬ 
portant topic, as a low oxygen mass fraction could be indicative 
of a small cross section of the poorly measured C 12 (a, y)0 16 nu¬ 
clear reaction, since it h as imp lications for the evolution of white 
dwarf progenitors (ISalaris et al.l[l997 i. Moreover, since the pre¬ 
cise value of the C 12 (a, y)0 16 cross section is s till the subject 
of an active debate - see, for instance, lAvila et aD (12015b . and 
references therein - any additional piece of evidence helping in 
constraining it is valuable. 

Furthermore, recent theoretical models predict that at moder¬ 
ately low luminosities a phase of stable hydrogen nuclear burn¬ 


ing in a shell ensues in the atmospheres of white dwarfs descend- 
ing f rom low metalli city progenitors dMiller Bertolami et aD 
12013b . More recently, Althaus et al.l (f2015l ) have expanded the 
calculations of iMiller Bertolami et al.l (120 1 3h and have shown 
that for progenitor metallicities between 0.00003 and 0.001, and 
in the absence of carbon enrichment due to the occurrence of 
a third dredge-up episode, the resulting hydrogen envelope of 
low-mass white dwarfs is thick enough to make stable hydrogen 
burning the most important energy source even at low luminosi¬ 
ties. Although in some cases this may not be a relevant source of 
energy, it may have non-negligible effects when computing the 
ages of very old clusters of low metallicity. 

In this paper we analyze the white dwarf population of the 
old, metal-poor globular cluster NGC 6397. We use an up- 
to-date population synthesis code based on Monte Carlo tech¬ 
niques, that incorporates the most recent and reliable cooling se¬ 
quences, a state-of-the-art description of the main properties of 
NGC 6397, as well as an accurate modeling of the observational 
biases. Our paper is organized as follows. In Sects. [2] and [3] we 
describe, respectively, our population synthesis code and the ob¬ 
served sample. Sect. [4] is devoted to analyze the results of our 
simulations. Specifically, we study the role of residual hydrogen 
burning, the effects of the slope of the adopted initial mass func¬ 
tion, the possible effects of mass segregation on the white dwarf 
population, and the fraction of white dwarf with hydrogen poor 
atmospheres. All this is done by performing a detailed statistical 
analysis of the luminosity function, of the color distribution and 
of the color-magnitude diagram of the white dwarf population of 
NGC 6397. Finally, in Sect. [5]our main results are summarized 
and our conclusions are drawn. 


2. The population synthesis code 

Our synthetic population code is based on Monte Carlo 
techni ques, and has been extensively used for st udying the 
disk (IGa rcfa-Berro et al.l 1 19991: Tbrres et al.l 1200 lb and halo 
(iTorres et al.ll2002 : IGarcfa-Berro et al.ll2004b populations of sin¬ 
gle white dwarfs. An improved version of this population syn¬ 
thesis code has been recently employed to model the prop¬ 
erties of the disk pop ulation of white dwarf plus main se¬ 
quence binary stars (ICamacho et al.l [2014b . Finally, and most 
relevant for this study, our code has also been successfully 
used to model Galactic open c lusters, such as NGC 6791 
(iGarcfa-Berro et all 120101 201 lbb. and globular clusters, like 
47 Tuc (IGarcfa-Berro et al. 1201 4b Detailed descriptions of our 
Monte Carlo simulator can be found in these papers. Thus, in 
this section we only summarize its most salient features, and we 
refer the interested reader for more details to the previously men¬ 
tioned works. 

Synthetic main sequence stars are randomly drawn accord¬ 
ing to a Salpeter-like initial mass function. In this_work we used 
the so-called “universal” mass function of iKrounal (12001) . For 
the mass range relevant to our study this initial mass function 
is totally equivalent to a two-branch power law with exponent 
-a, with a = 1.3 for 0.08 < M/Mq < 0.5 and a = 2.30 for 
M/Mq > 0.5. However, for the sake of completeness and sim¬ 
plicity we also use a classical Salpeter-like ( Salpeterl 1 1955b ini¬ 
tial mass function in Sect. 14.31 This mass function is a power law 
with just one index, which we consider a free parameter to fit the 
observations. The selected range of masses at the zero age main 
sequence is that necessary to produce white dwarf progenitors 
of suitable masses for NGC 6397. In particular, a lower limit of 
M > 0.5 Mq guarantees that enough white dwarfs are produced 
for a broad range of cluster ages. Within this mass range both 
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Table 1.x 2 test of the luminosity function, color distribution and color-magnitude diagram for different ages, duration of the burst 
of star formation and binary fractions. We list the reduced and normalized value of x 2 , that is the value of xl over its minimum 
value. 


x l test 


X F814W,v min 

Tf 606W-F8 14W.v 

/y 2 . 


x'Uxlin 







Tc = 

11.0 (Gyr) 




At (Gyr) 

1.0 

2.0 

3.0 

1.0 

2.0 

3.0 

1.0 

2.0 

3.0 

/bin = 

0.00 

9.42 

11.31 

12.68 

2.69 

3.60 

4.67 

4.17 

5.33 

5.88 

/bin = 

0.02 

9.72 

11.48 

12.79 

2.62 

3.54 

4.47 

4.18 

5.29 

5.90 

/bin = 

0.04 

9.71 

11.38 

12.86 

2.65 

3.56 

4.63 

4.16 

5.31 

5.97 

/bin = 

0.06 

9.62 

11.53 

12.84 

2.58 

3.54 

4.49 

4.16 

5.36 

5.82 






Tc = 

12.0 (Gyr) 




At (Gyr) 

1.0 

2.0 

3.0 

1.0 

2.0 

3.0 

1.0 

2.0 

3.0 

/bin = 

0.00 

4.43 

5.92 

7.33 

1.24 

1.68 

2.40 

1.71 

2.50 

3.28 

/bin = 

0.02 

4.49 

5.95 

7.54 

1.28 

1.62 

2.37 

1.78 

2.44 

3.23 

/bin = 

0.04 

4.44 

6.02 

7.52 

1.26 

1.67 

2.36 

1.72 

2.41 

3.20 

/bin = 

0.06 

4.38 

6.07 

7.39 

1.30 

1.63 

2.31 

1.71 

2.45 

3.13 






Tc = 

13.0 (Gyr) 




At (Gyr) 

1.0 

2.0 

3.0 

1.0 

2.0 

3.0 

1.0 

2.0 

3.0 

/bin = 

0.00 

1.95 

2.11 

2.86 

1.24 

1.08 

1.16 

1.43 

1.17 

1.48 

/bin = 

0.02 

2.05 

2.14 

2.92 

1.34 

1.07 

1.16 

1.42 

1.10 

1.45 

/bin = 

0.04 

1.90 

2.09 

2.94 

1.32 

1.08 

1.15 

1.38 

1.12 

1.41 

/bin = 

0.06 

1.87 

2.07 

2.97 

1.34 

1.11 

1.17 

1.37 

1.13 

1.40 






Tc = 

14.0 (Gyr) 




At (Gyr) 

1.0 

2.0 

3.0 

1.0 

2.0 

3.0 

1.0 

2.0 

3.0 

/bin = 

0.00 

4.00 

2.04 

1.18 

2.40 

1.75 

1.29 

2.59 

1.93 

1.34 

/bin = 

0.02 

3.81 

2.10 

1.11 

2.52 

1.82 

1.28 

2.60 

1.84 

1.39 

/bin = 

0.04 

3.72 

1.90 

1.08 

2.46 

1.78 

1.30 

2.61 

1.81 

1.28 

/bin = 

0.06 

3.61 

1.82 

1.07 

2.52 

1.79 

1.33 

2.54 

1.86 

1.38 


initial mass functions are totally equivalent when the standard 
value of o' = 2.35 is adopted. 

We adopt an age of the cluster, T c , which falls within an in¬ 
terval going from 11.0 to 14.5 Gyr. This is done to account for 
the widest relevant range of age estimates of NGC 6397. We also 
employ a star formation rate consisting of a burst of duration At. 
The width of the burst covers a sufficiently large range of values, 
from 0.0 to 4.0 Gyr, to ensure that all possibilities are explored. 
Finally, we adopted a maximum fraction of unresolved double¬ 
degenerate binaries /bin < 8%. This upp er limit i s consistent 
with the expected values reported by lHansen et al.l (120071) and 
iDavis et ail d2008h . The masses of the individual components of 
the binary system were randomly drawn according to a flat dis¬ 
tribution. However, the results are not sensitive to the choice of 
the mass distribution. We also considered that half of these bi¬ 
nary systems are composed by a helium-core white dwarf and a 
carbon-oxygen one, whereas for the rest of the systems the two 
components are typical carbon-oxygen white dwarfs. For white 
dwarfs with helium cores we adopted a mean mass of 0.4 Mq 
with a Gaussian deviation of 0.15 Mq 

Once we know which stars had time to evolve to white 
dwarfs we interpolate their photometric propert ies u sing a set 
of theoretical cooling sequences dMiller Bertolami et alil2013t 
[Althaus et al.ll2015 ) for white dwarfs with hydrogen-rich atmo¬ 
spheres, in accordance with the observed spectroscopic classifi¬ 
cation of the white dwarfs of the cluster. This set of evolution¬ 
ary sequences consists of a grid of cooling tracks for carbon- 
oxygen white dwarfs derived from evolutionary calculations of 
progenitors with metallicities ranging from Z = 1.0 x 10 4 to 
Z = 5.0 x 10 4 . We note here that these evolutionary sequences 
were evolved self-consistently from the zero age main sequence, 
through the giant phase, the thermally pulsing asymptotic gi¬ 
ant branch (AGB) and mass-loss phases, and ultimately to the 
white dwarf stage. We also emphasize that these evolutionary 
sequences provide us with a self-consistent initial-to-final mass 
relationship. For more massive oxygen-neon white dwarfs we 


employed the coolin g sequences of lAlthaus et al.l (120051) and 
lAlthaus et alJ (120071) . Fo r helium-core white dwarfs we used 
the cooling sequences of ISerenelli et alJ (1200 ll) , while for non- 
DA w hite dwarfs we us ed the cooling tracks of iBergeron et al.l 
(i201 ll) . and the colors of lSirianni et al.l (120051) . An important de¬ 
tail here is that to keep consistency, the masses of the result¬ 
ing synthetic white dwarfs were interpolated using the initial-to- 
final mass relationships obtained from the previously mentioned 
set of full evolutionary sequences for metal poor progenitors. 
Nevertheless, for the sake of completeness, we also ran some 
simulations for which the semi-e mpirica l i nitial -to-final mass 
relationship of field white dwarfs dCatalan et al.l l2008al lbh was 
used, obtaining essentially the same results. 

Photometric errors were assigned randomly according to 
the observed distribution. Specifically, for each synthetic white 
dwarf the photometric errors are drawn within a hyperboli- 
cally increasing band limited by cr\ - 0.2 0 «F 8 i 4 w — 31.0)~ 2 and 
cr u - 0.06 — 1.7 (ffiFRi 4 w - 31.0) -2 , which fits well the observa¬ 
tions of lHansen et al.! (l200 7i) for the /bfsmw magnitude. Similar 
expressions were employed for the rest of the magnitudes. 

3. The observed data set 

The observed data to which we compare ou r popul ation syn¬ 
thesis models was obtained from lHansen et al.l 1 20071) . and con¬ 
sists of a series of images taken with the Advanced Camera for 
Surveys (ACS) field located 5' southeast of the cluster core. Zero 
points for the instrumental magnitudes are 32.414 for wip 8 i 4 w 
and 33.321 for //ik, 06 w■ Besides, an absolute sharp parameter 
smaller than 0.02 and round parameter smaller than 0.02 was 
used to remove extended sources. In Fig. Q] we show the raw 
color-magnitude diagram of the sample. Those objects with 
bif814w < 2.3(mp606w - '«F8i4w) + 22.5 are far below the main- 
sequence and can safely be considered as candidate white dwarfs 
(black dots in the top-right panel ofFig.|T}. For magnitude 27 the 
completeness of the sample is 80%, and for magnitudes larger 
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Figure 1. White dwarf luminosity function, color-magnitude di¬ 
agram and color distribution of NGC 6397, for our best-fit 
model. Grey dots represent observed main sequence stars, black 
dots correspond to white dwarfs, while red points denote the re¬ 
sults of our simulations. The green squares represent the regions 
of the color-magnitude diagram for which we peformed ax 2 test, 
while the blue thin lines correspond to the cuts adopted to com¬ 
pute the distributions. The red curves correspond to the simu¬ 
lated distributions, while the blue ones are the observed distri¬ 
butions computed using our cuts. See the online edition of the 
journal for a color version of this figure, and the main text for 
additional details. 


than this value steadily decreases, reaching 50% at magnitude 
28. The black line in the top-left panel of this figure represents 
this limit. 

A third epoch set of observations to get deeper proper 
motions, which are especially useful for discriminating white 
dwarfs with faint magnitudes, has not been published yet. 
Consequently, and in order to avoid contamination from galax¬ 
ies, we apply a set of color cuts (represented by the blue lines 
in the color-magnitude diagram of Fig. [TJ. Those objects below 
the blue line of Fig.|T|are not considered in our analysis. Finally, 
we build a “clean” luminosity function and a color distribution 
function - left and bottom panels of Fig. Q] respectively - taking 
into account these additional color cuts (blue lines). 


4. Results 

4. 1. The color-magnitude diagram 

We first discuss the overall shape of the color-magnitude dia¬ 
gram, and the corresponding distributions of magnitudes and 
colors. All this information is displayed in Fig. [I] for a model 
in which we adopt an age T c = 12.8 Gyr, a burst duration 
At - 1.0 Gyr, and a fraction of binaries /bin = 4%, which we 
consider our reference model. The top-right panel of this figure 
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Figure 2. Probability density distribution for different pairs of 
free parameters (gray scale). Marked with a cross are the best- 
fit values, while the red lines represents the 68%, 90% and 95% 
confidence levels. The white cross corresponds to our reference 
model. 


shows the observed stars, as well as the synthetic white dwarfs. 
The results of our synthetic population calculations are repre¬ 
sented using red dots in the color-magnitude diagram, whereas 
the corresponding magnitude and color distributions are shown 
as red lines in Fig. [T] To obtain the color-magnitude diagram 
shown in this figure we fitted the position of the bright branch 
of the degenerate cooling sequence, obtaining an apparent dis¬ 
tance modulus (m - M )fsww = 12.42^^ and a color excess 
/AF606 W - F814W) = 0 2 2 ± 0.02. These values agree with 
those of lRicher et al.1 (2013) . and references therein. 

As can be seen, the agreement between our reference model 
and the observations is excellent. However, in our modeling 
we went one step beyond and performed a x 2 analysis using 
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the following strategy. We computed independent x 2 tests for 
the magnitude (Tf 814w ) and color CTf 606 w-f 814 w) distributions. 
Additionally, we calculated the number of white dwarfs inside 
each of the green boxes in the color-magnitude diagram of 
mm- which are the same regions of this diagram used by 
lHansen et al.l d2013l) to compare observations and simulations - 
and we performed an additional^- 2 test, xl>■ Then, we obtained 
for each of the three independent x 2 tests the reduced value of 
X 2 , xl = X 1 ! v ' where v is the number of degrees of freedom, 
v = N - p - 1, being N the number of bins and p the number of 
free parameters. Then the reduced;^ 2 values were normalized to 
the minimum value for each of the tests and added in quadrature. 
Employing this procedure only one estimator must be minimized 
to find the model that best fits all distributions. Clearly, this strat¬ 
egy is equivalent to a maximum likelihood method, hence the 
results of both approaches are consistently similar. For illustra¬ 
tive purposes, in Table[j]we list a representative set of values for 
the different tests and free parameters studied here. It is worth 
mentioning here that although we have used a x 2 test, our final 
aim is to estimate the values of the free parameters that best fit 
the observed data, rather than obtaining an absolute probability 
of agreement of our models with the observed data. 

In Fig.[2]we show using a gray scale the probability distribu¬ 
tion for different pairs of the studied free parameters - the cluster 
age, T c , the burst duration. At, and the binary fraction, /bin- Also 
shown, using red lines, are the curves enclosing the regions of 
68%, 90% and 95% confidence level. Finally, our best fit model 
is marked with a red cross and our reference model with a white 
one. Using the procedure previously detailed we found that the 
model that best fits the observed data has the following prop¬ 
erties. The age of the cluster is T c — 13.9 Gyr, the duration of 
the burst of star formation is At = 3.3 Gyr and the binary frac¬ 
tion of the cluster is /bin = 4%. One point of concern regarding 
our best-fit model is its very large age spread, 3.3 Gyr. This is 
un usua l in mos t globular cluster but on two, u> Centauri - see 
IVillanova et alj (j2014 f for_a recent discussion about this issue - 
and M54 (fslegel et al]|2007t) . Nevertheless, we remark that the 
previous theoretical efforts to model the white dwarf luminosity 
function of NGC 6397 could only achieve satisfactory fits when 
a constant star for mati on rate during the entire life of the cluster 
was employed dWinget et al.ll2009 ). However, we judge that the 
burst duration obtained using this procedure is totally unrealis¬ 
tic, and that our reference model, for which we adopted a typical 
burst duration of 1.0 Gyr is more realistic. We elaborate on the 
reasons for this below. 

Inspection of the top and bottom panels of Fig.[2]reveals that 
there is no correlation between the adopted fraction of binaries 
and the duration of the burst or the cluster age. Accordingly, 
/bin can be considered as an independent free parameter. Thus, 
its value has no major effects on the final results. In contrast, 
the middle panel of Fig. [2] unveils a clear unphysical correla¬ 
tion between the duration of the burst and the age of the cluster. 
Actually, ages shorter than that of our best-fit model are possible 
provided that smaller values of the burst duration are adopted, 
and the reverse is also true. In fact, ages as short as 12.8 Gyr are 
compatible, at the 95% confidence level, with the observations if 
the initial burst of star formation lasted for 1.0 Gyr. Even a clus¬ 
ter age of 12.4 Gyr is compatible at the 90% confidence level 
with the observed data if the duration of the burst of star forma¬ 
tion is 0.1 Gyr. Moreover, if At = 0.0 Gyr is adopted the age 
of the cluster turns out_to be 12.3 4 Gyr, in agreement with 
other ages estimates (lHansen et al.ll2013 f). All this is due to the 
very flat maximum of the probability distributions displayed in 

Fig. E] 


Since this is an important issue we explored which is the ori¬ 
gin of this unphysical correlation between the age of the cluster 
and the duration of the episode of star formation computed using 
the available white dwarf luminosity function. We found that the 
main reason for this is that the observed stars with mp 8 i 4 w > 26 
have large photometric errors. Thus, the precise position of the 
maximum of the observed white dwarf luminosity function can¬ 
not be determined with good accuracy. Hence, age determina¬ 
tions relying on the observed luminosity function are heavily in¬ 
fluenced by photometric errors. Additionally, such an extended 
episode of star formation would influence the morphology of the 
color-magnitude diagram of main sequence stars of NGC 6397. 
Nowadays with the high-accuracy phot ometr ic data collected 
with the Hubble Space Telescope (iMilone et al.ll2012h this age 
spread should be detectable, and this is not the case. In conclu¬ 
sion, we judge that the adopted photometric dataset is not yet ap¬ 
propriate for this kind of analysis, and that there are solid reasons 
to suspect that the duration of the initial burst of star formation 
is shorter than that obtained from our analysis. Consequently, 
from now on we stick to our reference model, in which a burst 
duration of 1.0 Gyr is adopted. 

Finally, we checked the sensitivity of our results to the dis¬ 
tance modulus and the reddening. Specifically, we varied the red¬ 
dening by +0.02, the nominal errors in its determination, and 
we found that the derived age for the cluster when our reference 
model is adopted varies by +0.02 Gyr. We repeated the same pro¬ 
cedure for the value of the distance modulus and we found that in 
this case the age of the cluster varies by ±0.3 Gyr. Furthermore, 
fixing the duration of the initial burst of star formation we ob¬ 
tained a conservative estimate of the error associated to the de¬ 
rived age of the cluster, and also to the fraction of binaries. We 
obtained <tt c =^'75 Gyr, and cr/- BIN As a final analysis 

we also studied the robustness of our results to the complete¬ 
ness of the observed sample. In particular, we discarded 20% of 
synthetic white dwarfs with magnitudes smaller than 26, while 
for magnitudes between 26 and 28 we assumed that the com¬ 
pleteness decreases linearly until 50%, and removed the corre¬ 
sponding synthetic white dwarfs. We then computed the white 
dwarf luminosity function and color distribution of the resulting 
sample, and we found that the results were essentially the same. 
We thus conclude that the age derived for our reference model is 
robust. 

4.2. The role of residual nuclear burning 

As mentioned earlier, residual nuclear burning in white dwarf 
atmospheres can be an important source of energy, specially 
for white dwarfs descending from progenitors of low metal- 
licity. This is so because, in the absence of third dredge-up, 
low-metalli city progenitors depart fro m the AGB with thicker 
envelopes (llben & MacDonald! 11986b . thus resulting in white 
dwarfs with thicker hydrogen envelopes. Consequently, residual 
hydrogen burning is expected to become relevant in white dwarfs 
with low-metallicity progenitors. In particular, iRenedo et al.l 
(2010) found that in white dwarfs resulting from progenitors 
with sub-solar metallicity, residual hydrogen burning may ac¬ 
count about 30% of the luminosity when cooling has proceeded 
down to luminosities ranging from L ~ 10 ~ 2 Lq to 10 3 Lq. 
This contribution is even more relevant when very low metal¬ 
licity progenitors (Z = 0.0001) and white dwarf masses smaller 
than ~ 0.6 Mq are considered. In these cases nuclear reactions 
are the main contributor to the stellar lumin osity for luminosi¬ 
ties as low as log (L/Lq) - -3.2 dMiller Bertolami et al.ll2013h . 
For instance, consider an otherwise typical white dwarf with 
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Figure 3. Same as Fig.Qjbut for the model with no residual hy¬ 
drogen burning, see Sect. 14.21 


A7 W n = 0.6 Mq, descending from a progenitor with the metal- 
licity of NGC 6397. For this white dwarf it takes about 1.04 Gyr 
to cool down to log (L/Lq) - -3.2 when no nuclear burning is 
considered, while when residual nuclear burning is considered, 
the cooling age is ~ 1.49 Gyr. Although for fainter luminosi¬ 
ties these differences become smaller, residual hydrogen burn¬ 
ing might still play a significant role in shaping the white dwarf 
luminosity function at moderately low luminosities. Given the 
distance modulus of NGC 6397 these luminosities correspond 
to apparent magnitudes of about wfsmw ~ 24 ~ 26, in the hot 
branch of the white dwarf luminosity function. 

To assess the possible effect of residual nuclear burning we 
used a second set of cooling sequences where this source of en¬ 
ergy was disregarded. As before, we built a grid of sequences 
for different masses and ages and interpolated for the precise 
metallicity of NGC 6397, keeping constant the duration of the 
initial burst of star formation, for which we adopted, as previ¬ 
ously discussed At - 1.0 Gyr. Afterwards, we ran our Monte 
Carlo simulator for the full range of the free parameters and ob¬ 
tained the combined^ 2 for each simulation, employing the pro¬ 
cedure explained in Sect. 14.11 The white dwarf luminosity func¬ 
tion, the color-magnitude diagram and the color distribution for 
the best-fit model is shown in Fig. [3 Although the general agree¬ 
ment between the simulated and observed samples for the differ¬ 
ent distributions of Fig. [3] is again quite good, the fit of some 
features of the observed distributions is worse when compared 
with that obtained when residual nuclear burning is taken into 
account. For instance, as can be seen in Figs. Q] and [3]the small, 
secondary peak of the observed white dwarf luminosity func¬ 
tion at /«F 8 i 4 w ~ 26.75 is best fitted when the sequences includ¬ 
ing residual hydrogen burning are considered. Moreover, when 
no residual nuclar burning is considered an excess of simulated 
white dwarfs redder than »Jf606w-f814w ~ 1-2 is quite apparent 
in the bottom panel of Fig. [3] According to these findings, for the 
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Figure 4. Same as Fig. |T| but for the model in which the slope 
of the initial mass function is allowed to vary (see Sect. 14.31 for 
details). 


rest the calculations presented in this work we will only employ 
evolutionary sequences in which residual hydrogen burning in 
white dwarf atmospheres is considered. 

4.3. The initial mass function 

The initial mass function is a key ingredient to model the ob¬ 
served properties of globular clusters. The most commonly used 
initial mass function consists of a power law with a character¬ 
istic power index a, oc M~ a , where M is the mass at the 

zero age main sequence. It is generally accepted that the value 
of the power law index should not differ much from a = 2.35 
t Salpeteri l l955l l. However, it has been recently suggested that 
the init ial mass function of NGC 6397 could have a flatter slope 
- see ( [Richer et al.|[2008l ) and references therein. In this section 
we explore such a possibility. 

As we did previously we consider the model described in 
Sect. 14. II but now, we consider a to be a free parameter of the 
fit, we allow it to vary between 1.35 and 3.35, and as we did 
before we compute the corresponding combined value of x 2 . In 
Fig- G] we show the result of this procedure. The best-fit model 
corresponds to a — 1.95. The rest of the values of the free pa¬ 
rameters for this best-fit model are similar to that obtained when 
a standard Salpeter-like initial mass function with a - 2.35 is 
adopted - see Sect. 14.11 Nevertheless, including a as a new free 
parameter does not substantially improve the fit. Therefore, we 
conclude that, unless a more restrictive observational data set is 
employed, a standard value of a describes adequately the global 
observed properties of the white dwarf population of NGC 6397. 
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Figure 5. Same as Fig. [7] but for a model in which we adopt 
a fraction of non-DA white dwarfs of 20% (see Sect. 14.41 for 
details). 


4.4. The fraction of non-DA white dwarfs 

Although observations indicate that the fraction of non-DA 
white dwarfs in NGC 6397 is small, of the order of / n 0 n-DA ~ 4% 
( Hansen et al . 2007 fc IStrickler et al.ll2009 : ). it is worth exploring 
which are its effects on the color-magnitude diagram, and on the 
corresponding white dwarf luminosity function and color distri¬ 
bution. To this end, in Fig. [5] we show the results when an un¬ 
realistic fraction of 20% of non-DA white dwarfs is adopted. In 
this color-magnitude diagram non-DA white dwarfs are shown 
as blue dots. When fractions smaller than 20% are used very few 
non-DA white dwarfs are found. Thus, this figure should be re¬ 
garded as an extreme case to illustrate the effects of varying the 
fraction of non-DA white dwarfs. 

As previously done in the preceding sections we consider the 
model described in Sect. l4.1l but we now consider f non -DA to be 
a free parameter of the fit, and we compute the corresponding 
combined value of x 1 , keeping fixed now the fraction of unre¬ 
solved binaries to 0.04, the best fit value obtained previously. 
The value of the fraction of non-DA white dwarfs that best fits 
the observed data is / non -DA = 0. Hence, our analysis concurs 
with observations that NGC 6397 has a negligible fraction of 
non-DA white dwarfs. However, we note that small fractions of 
non-DA white dwarfs are allowed by the present set of observed 
data at the 95% confidence level. 

4.5. Mass segregation in NGC 6397 

Another important issue worth investigating is the dynamical 
evolution of NGC 6397. The dynamical evolution of globular 
clusters is primarily governed by the gravitational attraction be¬ 
tween individual stars, although the interplay of the evolution 
of their individual stars and binary encounters also plays an im¬ 



Figure6. Top panel: Clean white dwarf luminosity function of 
NGC 6397 (black line) compared with the synthetic luminos¬ 
ity function for the entire sample of white dwarfs (grey line), 
and that obtained after eliminating massive and low-mass white 
dwarfs (red line). Bottom panel: mass distribution as a function 
of the magnitude for our simulated sample (black dots). Also 
displayed in this panel, as horizontal red lines, are the mass cuts. 

portant role. The high stellar densities in the cores of globular 
clusters favor binary encounters. In these dynamical interactions 
the more massive star usually transfers kinetic energy to the less 
massive one. Eventually, high-mass stars clump in the core of 
the cluster, while low-mass ones move to the outskirts of it, or 
even escape the cluster. 

Given that the observed sample of white dwarfs employed 
in our study was selected from the ACS field, which is lo¬ 
cated 5' SE of the cluster center, it is reasonable to expect that 
mass segregation may play a role. Alscr the evaporation time 
of NGC 6397 is ~ 3 Gyr ( I lev 1 et al. 2012) . implying that an 
important fraction of stars, in particular those with low masses, 
could have escaped from the cluster. However, the existence of 
mass segregation in NGC 6397 is s till controversial - see, for 
instance, the recent paper of iMartinazzi et all (120141) . and refer¬ 
ences therein - and a study employing white dwarf stars would 
help in settling this issue. Although an exhaustive asssessment of 
the dynamical evolution of NGC 6397 using the properties of the 
white dwarf population is beyond the scope of the present paper, 
a simple approach is worth pursuing, since a comparison of the 
properties of the simulated population with the observed sample 
could shed some light on the existence of mass segregation. 

To test how mass segregation affects the white dwarf popu¬ 
lation of NGC 6397 we proceed as follows. The bottom panel of 
Fig- [6]displays the distribution of white dwarf masses as a func¬ 
tion of the magnitude for our synthetic population. Black points 
denote typical carbon-oxygen white dwarfs, whereas the grey 
points represent white dwarfs with masses larger than a vari¬ 
able threshold M u (upper red line) and smaller than a given ad- 
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justable value M\ (lower red line). We also show the observed 
white dwarf luminosity function, which is shown as a black 
line in the top panel of Fig. [6] while the grey line is used to 
represent the luminosity function of the entire synthetic sam¬ 
ple. Finally, the white dwarf luminosity function for stars with 
masses M\ < M < M u is represented using a red line. As can 
be seen in Fig. [6] the bright portion (mp 8 t 4 w - 23 ~ 26) of the 
luminosity function is dominated by low-mass white dwarfs. By 
varying the lower mass cut we fitted the shape of the hot branch 
of the luminosity function, and obtained that the value of M\ that 
best fits the luminosity function is 0.535 Mq. On the other hand, 
massive white dwarfs dominate the region beyond the maximum 
of the white dwarf luminosity function. By adjusting the value of 
the upper cut in masses the height of the maximum of the white 
dwarf luminosity function can be fitted. The upper mass limit 
obtained in this way is 0.905 Mq. These two values are those 
shown in Fig.[6]as horizontal red lines. Note that the white dwarf 
luminosity function obtained after applying these two mass cuts 
(the red line in the top panel of Fig. [6]) fits extremely well the ob¬ 
served luminosity function. In fact, for the particular model ana¬ 
lyzed in Fig. [6j which corresponds to a cluster age of 12.8 Gyr, a 
burst duration of 1.0 Gyr and a binary fraction of 10%, the value 
of the reduced^- 2 is 1.20 when the mass cuts are applied. This 
number has to be compared with that obtained when no mass 
cuts are employed, which is x 1 = 1 -67, indicating a much bet¬ 
ter fit. Obviously, these cuts reduce size of the synthetic white 
dwarf sample. For the best-fit mass cuts this reduction amounts 
to a 10%. Hence, we conclude that mass segregation plays a role 
in the fine tunning of the luminosity function. 


5. Conclusions 

In this paper we studied the white dwarf population of 
NGC 6397. We have done that employing an up-to-date Monte 
Carlo population synthesis code, which incorporates the most 
advanced descriptions of the observational biases and selection 
procedures, and state-of-the-art cooling sequences, for the ap¬ 
propriate metallicity of this cluster. We have compared the re¬ 
sults of our simulations with the most recent observational sam¬ 
ple of the white dwarf population of this cluster, which was ob¬ 
tained using deep Hubble Space Telescope observations. In par¬ 
ticular, our simulations were compared with the general appear¬ 
ance of the color-magnitude diagram, with the color distribution 
and with the white dwarf luminosity function in a quantitative 
manner, obtaining probability distributions for each pair of the 
most important free parameters of the fit. These are the age of 
the cluster, the duration of the initial burst of star formation, 
and the fraction of unresolved binaries. Also, other key ingredi¬ 
ents in modelling the white dwarf population of NGC 6397 were 
studied. These include the adopted reddening and distance of 
the cluster, the role of the completeness of the observed sample, 
the fraction of unresolved binaries, the initial-to-final mass rela¬ 
tionship, the adopted white dwarf cooling sequences, the initial 
mass function, the fraction of non-DA white dwarfs, and the im¬ 
pact of mass segregation. Our results are in very good agreement 
with the observed data, and are largely independent of model as¬ 
sumptions. However, the probability distributions show broad, 
flat maxima that prevent to simultaneously obtain accurate de¬ 
terminations of the age of the cluster, and of the duration of its 
initial burst of star formation. The origin of these flat maxima is 
the still sizable size of photometric errors of the current observed 
data set. 

The age of NGC 6397 derived using the white dwarf cool¬ 
ing sequence is 12.8^ ^ Gyr, when the duration of the initial 


Table 2. Age determinations for NGC 6397. 


Tc (Gyr) 

12.0 ± 0.8 

13.4 ±0.8 

13.9 ± 1.1 

11.47 ±0.47 

12 0 +0 ' 5 
- 1.0 

13.50 ±0.50 
13.00 ±0.25 


12.8 


+0.5 

-0.75 


Reference 

Anthony-Twarog & Twarog (2000) 

i Chabo^eretah (200T) 

^rattonjt^aL (200T) 

^fanser^Uif (2007) 

Wingete^ak (2009) 

^ottei^Uih (2010) 

VandenBer g et al. (2017) 

This work 


burst of star formation is 1.0 Gyr. The uncertainty in the age of 
the cluster introduced by the uncertainty in the distance modu¬ 
lus is +0.3 Gyr, while that introduced by the reddening is neg¬ 
ligible when compared with the previous one. This age estimate 
agrees well with the age determinations obtained using the main 
sequence turn-off - see Table [2] Specifically, our age determi¬ 
nation is slightly shorter than the most recent d eterm inations ob¬ 
tained using main sequence stars. Speci ficallv.lVandenBerg et~akl 
(120 1 3h obtained 13.00 ± 0.25 Gyr, and iDotter et al. ( 201 Oh de- 
rived 13.50 ± 0.50 Gy r, some what older than previous es¬ 
timates. For instance, GrattpngtalJ (120031) obtained 13.9 ± 
1.1 Gyr and IChabover et ah (2001) derived 13.4 ± 0.8 Gyr. 
However, our age estimate is in line with the age derived by 
[Anthonv-Twarog & Twarori(20001) . whoobtainedl2.0±0.8Gyr. 
We also find that ou r age estimat e for NGC 6397 is significantly 
older than th o se of Hansen et alJ d2007l) . 11.47 + 0.47 Gyr, and 
IWinget et al.l (12009 ). 12.Oj^ Gyr, which were obtained using 
the position of the cut-off of the white dwarf luminosity func¬ 
tion. This discrepancy can be partially attributed to the fact that 
we used updated evolutionary tracks. This includes not only the 
white dwarf cooling sequences, but also the main sequence life¬ 
times of white dwarf progenitors, and the initial-final mass re¬ 
lationship, which in our calculations are computed for the ap¬ 
propriate metallicity of NGC 6397. However, the use of reliable 
evolutionary tracks is not the only reason for this, and we stress 
that other factors, like the way in which the synthetic population 
of white dwarfs is generated, taking into account all the known 
observational biases and selection procedures, are also important 
to obtain reliable ages. 

Another interesting result of our study is that we found that 
the fraction of binaries is close to zero. This is particularly re¬ 
assuring, as the percentage of binaries in our best-fit model 
(~ 4.0%) is very simil ar to tha t obta ined for main sequence stars 
- see tables 2 and 3 in lMilone et al.l d2012l) . 

Having estimated the age and the percentage of unresolved 
binaries of NGC 6397 we decided to study other important prop¬ 
erties of the model. To start with, we studied the role of hydro¬ 
gen residual burning in the atmospheres of white dwarfs with 
very metal poor progenitors. To this end, we computed a com¬ 
plete set of full evolutionary sequences in which residual nuclear 
burning was (artificially) not considered. The results of this nu¬ 
merical experiment indicate that although good fits to the ob¬ 
served properties of the white dwarf population of NGC 6397 
can also be obtained using this set of evolutionary tracks, the fit 
is better when the cooling sequences in which residual hydrogen 
burning is considered. This result sheds light on the importance 
of third dredge-up and extra mixing episodes in low-mass, low- 
metallicity AGB stars, for which theoretical and observational 
evidence js notconclusive about their occurrence. As shown in 
lAlthaus et al.l (1201 5l) . in the absence of third dredge-up episodes 
during the AGB phase, most of the evolution of white dwarfs re- 
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suiting from low-mass (M < 1.25 Mq ), low-metallicity progen¬ 
itors is dominated by stable hydrogen burning. In view of our 
result that low-mass white dwarfs in NGC 6397 are expected 
to sustain significant residual hydrogen burning, we conclude 
that their low-mass progenitors would not have suffered from 
the carbon-enrichment due to third dredge-up during their AGB 
evolution. 

Another possible concern could be the slope of the initial 
mass function. To assess this point we conducted a series of sim¬ 
ulations in which we varied the power-law index of the initial 
mass function, and we found that that although a shallower slope 
will fit slightly better the properties of the observed sample, the 
standard value of the Salpeter mass function fits equally well 
such properties, within a 95% confidence level. Also, we studied 
which is the impact of adopting a different fraction of non-DA 
white dwarfs and found that our best fit model corresponds with 
a negligible fraction of these stars, in agreement with observa¬ 
tions. Finally, we made a preliminary test to assess the effects of 
mass segregation on the white dwarf population of NGC 6397. 
By adjusting the bright and the dim portion of the luminosity 
function we obtained that white dwarf with masses larger than 
0.905 Mq or smaller than 0.535 Mq are absent in the observed 
sample, possibly as a consequence of mass segregation. When 
this is taken into account in the simulated populations the agree¬ 
ment between the observed and the theoretical data turns out to 
be excellent. 

In summary, we have conducted the most complete and up- 
to-date population synthesis study of the white dwarf population 
of the old, metal-poor, globular cluster NGC 6397. This study 
has allowed us to derive useful constraints on the characteristics 
of this cluster, and to obtain an estimate of its age. Nevertheless, 
we emphasize that a better observational data set, with smaller 
photometric errors for white dwarfs of magnitudes ^ 27 will 
help in obtaining even more reliable estimates of both the age of 
its star formation history. 
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